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Abstract 

We investigate two-point correlation functions of left-handed currents computed in 
quenched lattice QCD with the Neuberger-Dirac operator. We consider two lattice 
spacings a ~ 0.09,0.12 fm and two different lattice extents L ~ 1.5,2.0 fm; quark 
masses span both the p- and the e-regimes. We compare the results with the predic- 
tions of quenched chiral perturbation theory, with the purpose of testing to what extent 
the effective theory reproduces quenched QCD at low energy. In the p-regime we test 
volume and quark mass dependence of the pseudoscalar decay constant and mass; in 
the e-regime, we investigate volume and topology dependence of the correlators. While 
the leading order behaviour predicted by the effective theory is very well reproduced 
by the lattice data in the range of parameters that we explored, our numerical data are 
not precise enough to test next-to-leading order effects. 
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1 Introduction 



Lattice QCD underwent continuous advances in recent years, thanks to increasing com- 
puting resources and to very important theoretical and algorithmic improvements. One 
of the most important outcomes of these efforts is that unquenched computations arc 
now reaching ranges of volumes and quark masses where the QCD chiral dynamics can 
be tested and the matching with the chiral effective theory be performed. Several stud- 
ies for Nf = 2 have been recently undertaken, with Wilson fermions (plain and with 
0(a) improvement [1,2]) and Wilson twisted mass fermions [3-5]. For the case Nf = 3, 
beside the results published in 2004 obtained with staggered fermions [6] (see [7] for 
recent updates), new data are available for Domain Wall fermions [8] and for Wilson 
fermions [9]. Simulations with pseudoscalar masses Mp as low as 200 — 300 MeV are be- 
coming state-of-the-art. Therefore, it is highly important to understand to what extent 
the chiral effective theory can be applied in the region, in order to perform controlled 
chiral extrapolations. 

On the other hand, in the past years many studies have matched quenched lattice results 
with the chiral effective theory [10-14], using the non-trivial assumption that quenched 
QCD at low energy is described by quenched chiral effective theory. The main goal of 
this work is to verify this assumption in two different kinematic corners of the chiral 
regime of QCD. 

On a finite volume with linear extent L, apart from the conventional p-regime, where 
finite- volume effects are exponentially suppressed in (MpL), one can extract information 
from other kinematic regions, for instance the e-regime, where the chiral limit m — > 
is taken while keeping the pion wavelength much larger than L [15, 16]. Higher order 
corrections in the chiral theory for physical observables look very different in the two 
cases; the matching between QCD and the chiral effective theory gives rise to different 
systematic effects and allows to extract quantities in the chiral limit from different 
observables. On the lattice, several quenched simulations in the e- regime with Ginsparg- 
Wilson fermions have been performed [13,17-24]. In the quenched approximation it is 
now possible to collect results with Ginsparg- Wilson fermions at different volumes and 
lattice spacings with high statistics and reasonable computational effort. On the other 
hand, in full QCD the simulation of the Ginsparg- Wilson Dirac operator is still very 
expensive from the numerical point of view. First results for Nf = 2 are available both 
in the p- and e-regime [25-32] , although at the moment only for limited ranges of lattice 
spacings and volumes. For other regularisations, which break explicitly chiral symmetry 
at finite lattice spacing, exploring the e-regime is problematic. 

We will adopt the Neuberger solution of the Ginsparg- Wilson relation for the Dirac 
operator. In particular we will compute the left-handed current two-point correlation 
function both in the p- and in the e-regimes; in the p-regime we will investigate the 
quark mass and volume dependence of the pseudoscalar decay constant and mass, and 
compare it with the expectations from the chiral effective theory at NLO. In the e- 
regime, we will also study the topology dependence of the current correlators. We 
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shall extract the corresponding Low Energy Couplings (LECs) and compare the results 
obtained for the leading order constants E and F in the two regimes. Moreover, by 
considering different lattice spacings, we will check that lattice artifacts are small, as 
already observed in many computations using the Neuberger Dirac operator (see e.g 
[33]). This work complements and expands the study published in [13,33]. 
The paper is organised as follows: in Section [2] we recall the main results from the 
quenched chiral effective theory at NLO, in the e- and p-regime; Section [3] is devoted 
to describing the details of our numerical simulations; in Sections [U and [5] we match 
the lattice results with the chiral effective theory in the e- and p-regimes respectively; 
finally, in Section [6] we compare and discuss the results obtained for F and E in the two 
cases. 



2 Current correlator in the chiral effective theory 

At leading order, the Euclidean Lagrangian of the chiral effective theory is given by 
[34, 35] 

£ = El Tr [d^djj] - |Tr {eWtUM + MWe-"' N *} , (2.1) 

where U £ SU(Nf) contains the pseudoscalar degrees of freedom and M is the mass 
matrix. For simplicity we consider a mass matrix proportional to the identity, A4 = mi. 
F and S are the pseudoscalar decay constant and the quark condensate in the chiral 
limit, and 8 represents the vacuum angle. 

In a finite volume V = TL 3 with L S> 1/Aqcd ; one can distinguish different 
chiral regimes. If Mp is the pseudoscalar meson mass, approaching the chiral limit by 
keeping MpL 3> 1 defines the so-called p-regime. In this case the chiral effective theory 
looks essentially as in the infinite volume case: finite-volume effects are exponentially 
suppressed by factors exp(— MpL), while the mass effects are the dominant ones. The 
power counting in terms of the momentum p and quark mass m is given by 

m~p 2 , 1/L,l/T~p. (2.2) 

Alternatively, one can approach the chiral limit while keeping n = raEF < 1; in this case 
the Compton wavelength associated with the pseudo-Goldstone bosons is much larger 
than the linear extent L of the box, and volume effects are enhanced. This defines the 
e-regime, where the power-counting is reorganised such that [15, 16] 

m~e 4 , 1/L,l/T~e. (2.3) 

One of the most important effects of the reorganisation of the power counting is that, at 
a given order in the effective theory, fewer Low Energy Couplings (LECs) appear with 
respect to the p-expansion. The fact that the corresponding higher-order counterterms 
are kinematically suppressed may be convenient for the extraction of LECs by matching 
the effective theory to lattice QCD. 
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In this work we consider the left-handed current, which at leading order in the 
effective theory formalism corresponds to 

J£ = ^Tr (T a Ud,U^ , (2.4) 

where T a are the traceless generators of SU(iVf). In particular we are interested in the 
two-point correlation function 

C ab (t) = J d 3 x{J a {x)J b {0)) = Tr[T a T b ]C(t). (2.5) 

The chiral effective theory formalism can be extended to the quenched case; in partic- 
ular, two equivalent methods have been developed to cancel the fermion determinant, 
namely the supersymmetric formulation and the replica method [36-38]. An impor- 
tant feature of the quenched setup is that the flavour singlet does not decouple in this 
case; moreover, its mass parameter ffig/ (2N C ) is related to the topological susceptibility. 
In the following we summarise the known results for this current correlator from quenched 
chiral perturbation theory at NLO0, in the e- and p-regimes, with degenerate quark 
masses. Current correlators have been recently computed in the effective theory also 
for non-degenerate quark masses, in the full and partially quenched scenarios, in the 
case where all quarks are in the p- or in the e-regime and in the mixed case, where 
m v T,V < 1 for the valence quarks and m s TiV 3> 1 for the sea quarks [40,41]. 



2.1 



e-regime 



Topology plays a relevant role in the e-regime [42], such that observables must be consid- 
ered in sectors of fixed topological charge. In the quenched case, the current correlator 
in Eq. (12.5(> at NLO and fixed topology v is given by [43, 44] 



Cu(t) 



El 

2T 



1 + Jvy a M h l 



with 



and 



hi(r) 



1 

12 



ov(ju) = ix [I y (n)K v (n) + I v+1 (fi)K v -i(fj,)] + 



/i 



(2.6) 



(2.7) 



(2i 



where I u and K u are modified Bessel functions. The most notable fact is that in the 
NLO expression only the leading-order LECs S and F enter, as already anticipated. 
In our analysis we will compare both the time and topology dependence of the QCD 



1 We neglect terms proportional to which is the parameter associated to the kinetic term of the 
singlet field, since it is suppressed for this observable in the simultaneous expansion in momenta and 
1/JVc [39]. 
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correlators computed on the lattice with the expectations of the effective theory. A 
convenient way to study the topology dependence is to fix t = T/2; in the chiral limit 
one has 

^ov(/x)! M= o = M, (2.9) 

hence one expects C v {T/2) to depend linearly on the topological charge v. Moreover, 
one obtains the parameter-free prediction 

24L 3 [C vl (T/2) - C U2 (T/2)} | M=0 = H - N- (2.10) 

From this expression and from Eq. (|2.6f) it becomes clear that the sensitivity to topology 
is quite limited: in order for C u (T/2) to be significantly different from C u +/\ u (T/2) one 
needs the following condition on the relative error: 

AC ^ T / 2 ) « ^ 1 (2 n) 

Using the quenched value F ~ 100 MeV from [13], this implies that statistical errors 
much smaller than ~ [(14Au)T/L]% for L = 1.5 fm and ~ [(8Au)T/L)% for L = 2.0 
fm must be reached. Notice that NLO effects are larger for asymmetric boxes, since 
they are proportional to (T/L) 3 ; however, if T 3> L one enters in a different kinematic 
range, called 5-regime, which will not be discussed in this work. 
For /i<l, the leading /^-dependence in the NLO correction is given by 

""W^n^im u To! <2 ' 12) 

[ Ll-7-log(f)J n A + ... (i/ = 0), 

where 7 is the Euler-Mascheroni constant. For i^/Owe then expect a weak sensitivity 
to quark mass for the current correlator. 

Matching the left correlator computed in lattice QCD with the chiral effective 
theory allows to extract the low-energy constant F with control over NLO effects. For 
this particular correlator the chiral condensate S appears only at NLO; in particular the 
expression in Eq. (|2.8p represents the quenched chiral condensate at finite fi at leading 
order in the e-expansion [45,46]. At NLO, the condensate retains the same functional 
form of Eq. (|2.8p [45,47], with /u replaced by n e g = (mT, e gV) and 

S eff (F) = S [1 + w o H(0)] , (2.13) 



where 



Moreover, we have defined 



*W = V ^ 7^2 ( 2 - 14 ) 



2 
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where m^/(2N c ) is the flavor singlet mass parameter; as already anticipated, it is related 
to the topological susceptibility by the equatiorJl 



(u 2 ) _ m%F 2 



V 

In dimensional regularisation one obtains 



1 + 2ci + In 



where L = V 1 ^ 4 , 1/Lq is the ultraviolet subtraction point, and 



^2 + \(-7 + H**)). 



(2.16) 



(2.17) 



(2.18) 



@2 is a shape coefficient [48], which in the symmetric case T = L takes the value 

(3 2 = -0.020305. (2.19) 

The infrared "sickness" of quenched QCD is reflected here in the fact that £ e ff(V) 
diverges in the limit L — > oo. In the following we define operatively S e gf at a fixed V as 



2HE / n 

— tt <r v [p) 

mv \ fj, 



S c 2 ff (L), v > 0. 



(2.20) 



m=0 



2.2 p- regime 



In the p-regime, the NLO finite-volume prediction for the current correlator in the 
quenched case is given by [49] 



C { t) = h I UF^ C0Sh[{T/2 -' )Mi:] 



2sinh[TM^/2] 

The pseudoscalar decay constant in this case is volume-independent: 

M 2 



F P 



1 + 



r«5 



(2.21) 



(2.22) 



2(4vrF) 2 

where a, are the LECs associated with NLO operators in the (quenched) chiral La- 
grangian in the convention of [10], and 

2mS 



M 



F 2 



For the finite- volume pseudoscalar meson mass one obtains 

rV\ 2 _ A j-2 



{M V pY 



M P [l + w g 2 {M P ,V)} 



Mi 



M 1 



l + w H(M 2 ) 



M 2 



(4vrF) 



r(a 5 - 2a 8 ) 



(2.23) 

(2.24) 
(2.25) 



2 We adopt the normalisation conventions of [43,44]. 
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where H(M 2 ) is given, in dimensional regularisation, by 



H(M 2 



d d p 



1 



1 



(2ir) d (p 2 + M 2 ) 2 (An) 2 
The volume-dependent function g r reads [48] 



g r (M P ,V) 



1 



(4vr) s 



dX 



, M 2 



(2.26) 



(2.27) 



nez 4 



exp 



In our analysis we will investigate finite-volume effects by comparing lattice results 
obtained at different volumes V\ and V 2 . A convenient quantity to consider for this 
purpose is the ratio Mp 1 /Mp 2 , for which one obtains the NLO expression 

2 



Mp 1 
Mp 2 



92(M P ,V 1 ) - g 2 (Mp,V 2 ) 



(2.28) 



If wo is given as input, this is a parameter- free prediction from the chiral effective theory. 
By reabsorbing the divergences in Eq. (|2.26|) in the low-energy constant £ one obtains 



Ml 
2m 



F 2 



1 



w i M' 
log _ 



(4vr 



M 1 



(AnxFf 



(a 5 - 2a 8 ) 



(2.29) 



E(//) is related at NLO to the S e g at a given scale (-L e ff) defined in the e-regime in 
Eq. QZM) by 



^es(L c s) = £0u) 1 + ^0 /?2 + 



1 



(4vr) 2 



l + log^/, 2 )) 



(2.30) 



3 Numerical simulations 



For our numerical study we have adopted the Neuberger-Dirac operator D [50,51] and 
the Wilson gauge action, on a box with volume V = L 3 T and lattice spacing a. We 
have computed the two-point function 



C ab (t) = Y^(JS(x)J b (0)) = Tr[T a T b ]C(t), 

X 

with the left-handed current 
where P± = (1 ± 7s)/2 and 



J a = 4>T a l0 P-i>, 



1> 



l 



l 



1 + s' 



(3.1) 
(3.2) 

(3.3) 
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The parameter \s\ < 1 has been chosen equal to 0.4. For complete definitions and 
conventions related to the Neuberger-Dirac operator the reader can refer to [13]. An 
advantage of using left-handed currents is that zero-modes of the Dirac operator do not 
contribute to the corresponding correlator: at finite volume, no divergences are present 
in the m — > limit. 

We have considered two sets of lattices, one dedicated to simulations in the p-regime 
and one for the e-regime. The parameters of the two sets are reported in Tables [I] and [2j 
For the p-regime we have chosen two different lattice spacings and two spatial extents 
L ~ 1.5 fm (pi and p3) and L ~ 2 f m (p2). The temporal extent is chosen to be T = 2L 
(pl,p2) or T = 3L/2 (p3). For the lattice p3 we did not perform new simulations and 
we used instead the data already presented in [13] . 

For the e-regime, we have chosen three symmetric lattices (el, e2, e3) with the same 
parameters used in a previous work [24] for the computation of the quark condensate. 
The quark masses here have been chosen such that (mV)/(ZgrQ) is constant for the 
lattices (el, e2, e3), where Zg is the renormalisation constant of the RGI scalar density 
[33,52] and tq ~ 0.5 fm [53]. Moreover, we have considered two additional asymmetric 
lattices, e4 and e5, with T = 2L. For the lattice e4 we used the measurements collected 
in a previous project [54]. 

Following [13,55], we applied the low-mode averaging technique in order to reduce large 
fluctuations induced by low-modes wave functions. In Tables [Hand [21 N\ ow indicates the 
number of low-modes which have been extracted. The computation of the topological 
index, the low-lying eigenvalues and the inversion of the Neuberger Dirac operator have 
been performed using the techniques in [56] . For all the lattices, the low eigenvalues have 
been computed with a 5% precision. The values of the scale r$ in lattice units [57], Zg, 
and the renormalization factor for the local left-handed current [33] for our couplings (3 
are listed in Table El 



lat 


p V 


^cfg 


am 




Pi 


5.8458 12 3 x 24 


475 


0.02, 0.03, 0.04, 0.06 


7 


P 2 


5.8485 16 3 x 32 


197 


0.02, 0.03, 0.04, 0.06 


20 


p3 


6.0 16 3 x 24 


113 


0.025, 0.04,0.06,0.08,0.1 


8 



Table 1: Simulation parameters for the p-regime 



7 



lat 


p V 


NK 


am 




el 


5.8458 12 4 


177, 313,221, 
126,62,35 


0.001, 0.003, 0.008, 
0.012,0.016 


20 


e2 


5.8458 16 4 


49,69,82,72 
50,54,38,32 


0.000316, 0.000949,0.00253, 
0.00380,0.00506 


20 


e3 


6.0 16 4 


131,231,178, 
96,44,16 


0.000612, 0.00184,0.00490, 
0.00735, 0.00980 


20 


e4 


5.8485 16 3 x 32 


151, 130, 125, 101,87 
86,66,52, 29 


0.002,0.003 


20 


e5 


5.8458 12 3 x 24 


22, 55, 50, 56,25 
21,22 


0.003 
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Table 2: Simulation parameters for the e regime. N£ indicates the number of config- 
urations for topologies \v\ = 0,1,2..., except for lattice e4, where the lowest index is 
\v\ = 2. 



p 


r /a 


Zj 


Z s 


5.8458 
5.8485 
6.0 


4.026(23) 
4.048(23) 
5.37(3) 


1.710(5) 
1.706(5) 
1.553(2) 


1.28(6) 
1.28(6) 
1.05(5) 



Table 3: Values of the Sommer scale r$/a [57] and the renormalisation constants of the 
left-handed current and the scalar density [33] for the values of (3 used in our study. 



4 Matching lattice QCD with the chiral effective theory: e-regime 
4.1 LO matching 

At LO in the e-expansion, the left correlation function is expected to be independent 
on time, mass and topology; from Eq. (|2.6[) one reads 

Cu(t) = ^. (4.1) 

In Fig. [T] we plot the renormalised dimensionless quantity 2T 'Z"jC u (T '/2)rg as a function 
of T/ro, for \u\ = 2 . We do not observe a significant mass dependence of this quantity, 
hence we report the results only at a single value of \x (see figure caption). 
The results obtained for lattices el (black filled circles) and e3 (blue filled triangles) 
corresponding to T ~ 3ro are in good agreement between each other; this indicates that 
within our precision we are not sensitive to lattice artifacts. Their consistency with 
the data of lattice e2 (red filled squares), for which T ~ 4ro, is a verification of the 
1/2T scaling predicted at LO. In this plot we also show the data obtained with the 
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M=2 



- 

_ . el 
H e2 

o 

a e3 

□ e4 -L 
o e5 

J i I , I , I , I , I 

3 4 5 6 7 8 

T/r Q 

Figure l: The quantity 2TZjC u (T/2)r% plotted clS cl function of T/r$ for lattices el 
(black filled circles), e2 (red filled squares) and e3 (blue filled triangles). The plot 
reports the data at fixed topological charge \v\ = 2 and fixed fj, ~ 0.6, corresponding 
to the intermediate quark mass. The red empty squares correspond to the lattice e4 
(fi ~ 0.5, corresponding to the lightest mass), while the black empty circles correspond 
to the lattice e5 (fi ~ 1). 

asymmetric lattices e4 (for am = 0.002, corresponding to fi ~ 1) and e5 (/i ~ 0.5). In 
this case the data tend to depart from the results of the symmetric lattices, although 
not significantly within our statistical errors. From the chiral effective theory we expect 
NLO effects to be larger for asymmetric volumes, as already discussed in section [2j 

4.2 NLO matching: topology dependence 

At NLO, the e-expansion gives predictions on time, mass, and topology dependence of 
C u (t); in particular, one expects C u (T/2) at the chiral limit to depend linearly on the 
topological charge. 

In Fig. M we show the differences 2AL 2 Zj [C Vl (T/2) - C V2 {T/2)) for several choices of 
v\ , z/2 and compare the results with the parameter- free predictions of the chiral effective 
theory, Eq. (|2.10p . For the lattices e2 and e5 we excluded the sector \v\ = 1 because 
of large statistical uncertainties. As already noticed in the previous section, we don't 
observe significant quark mass dependence. The statistical errors associated with the 
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24L Zj [ C V| (T/2) -Cv 2 (T/2) ] 
. el 



. e2 

A e3 

□ e4 

n e5 



2-1 



3-2 



3-1 



4-2 



7-2- 



5-1 6-2 



4-1 5-2 



f 



Figure 2: The differences 24L 3 Z5 [C Vl (T/2) - C U2 {T/2)) computed in the e-regime. The 
lines represent the parameter- free prediction from the chiral effective theory (Eq. (|2.10p ). 
and the numbers denote the different combinations of — \fi\). The data from lattices 
el, e2 and e3 refer to [i ~ 0.6, corresponding to the intermediate quark mass. For the 
lattice e4, the lightest quark mass fj, ~ 0.5 is considered. 



sector v = are very large. Hence, also in the case where we would expect a stronger 
quark mass dependence as predicted by Eq. (|2.12|) . we are not sensitive to mass effects. 
Anyway we exclude the v = sector from this analysis. 

Numerical data tend to depart systematically from the theoretical expectations when 
the topological charge is increased. Although statistical uncertainties associated with 
the results are too large for precise quantitative statements, this may indicate that 
higher order corrections are significant. 

4.3 NLO matching: time dependence and determination of F. 

At fixed values of the topological charge and quark mass, the e-expansion predicts a 
parabolic time dependence of current correlators. Since we have already observed that 
our sensitivity is not high enough to test NLO dependence, we do not enforce the 
topology and mass dependence of the NLO corrections and leave it as a free parameter. 
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In particular we perform the 2-parameter fit 



(4.2) 



The results for B\, B2 at the different masses and topological charges are collected in 
the Appendix [A] in Tables El [3 Three examples of the fit are given in Fig. [3l for 
the intermediate mass and \v\ = 1. The fit ranges are t/a = 5 — 7, t/a = 5 — 11 
and t/a = 6 — 10 for the lattices el, e2 and e3 respectively. Asymmetric lattices are 
not considered in this analysis. We observe that the results for B\ are stable within 
errors with respect to the quark mass and topology. Moreover, each topological sector 
yields an independent determination and it is possible to reduce the associated error by 
performing averages between different sectors. By varying the fit range, no significant 
deviation is observed for B\. We determine F in lattice units by averaging the results 
obtained for B\ at the five quark masses, and then performing an additional average over 
different topological sectors. In particular, for the lattices el, e3 we consider \u\ = 1 — 3, 
while for the lattice e2 we averaged in the interval \v\ = 2 — 4. This choice is due to 
large fluctuations which affect the sector \v\ = 1; the same happens for the v = sector, 
for all lattices. For the dimensionless quantity Fro we obtain 



The agreement between (el) and (e2) indicates that finite-volume effects are below our 
statistical precision; their agreement with (e3) is a signal that also lattice artifacts are 
smaller than our errors. 

In the chiral limit, the effective theory predicts the coefficient B2 to be equal to 
the topological index \v\\ in Fig. 0]we report the results for B2 at the smallest quark 
mass, as a function of \v\, together with the theoretical expectation. The results for 
B2 are more sensitive to a change of the fit range with respect to B\\ in particular 
for the lattice el, enlarging the fit range gives differences of the order of 2-3 standard 
deviations. This is not surprising, since the latter has a relatively coarse lattice spacing 
and a relatively small volume, and the number of points available for the fit is hence 
limited. Here we consider only the results concerning the fit range t/a = 5 — 7, keeping 
in mind that additional systematic error on B2 for lattice el can be substantial. 
Like in the previous analysis at fixed t = T/2, the data show a tendency to depart from 
the prediction with increasing indicating that higher order effects may be important. 
In particular, corrections are expected to be severe when 3> \J {v 2 ). However, also 
in this case the statistical errors are very large and do not allow for a precise statement. 
In this fit, the stability of the coefficient B\ with respect to m and v can be interpreted 
as a signal that systematic uncertainties coming from higher orders are under control. 



Fr 
Fr 
Fr 



0.284(4), (el) 
0.278(6), (e2) 
0.280(5). (e3) 



(4.3) 
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Figure 3: Time dependence of the correlator ZjC u (t) for \u\ = 1 at the intermediate 
quark mass, for the lattices el, e2, e3. The curves represent the fit of Eq. (|4.2|) . with 
ranges t/a = 5 - 7 (el), t/a = 5- 11 (e2) and f/a = 6 - 10 (e3). 



As additional test we can constrain the NLO terrd_| I?2 = \v\ 

2 

^c,(t) = ^ + jj hl (f) ■ (44) 

The results obtained for B\ for different \v\ are shown in Fig. [5] for the smallest quark 
mass; compatible results are obtained for the heavier quark masses. In this case the 
importance of higher order terms in the chiral theory manifests itself in the fact that 
B\ is topology dependent . However, for low \v\ we observe a constant behaviour 
within statistical errors: moreover, data at low \u\ are fully compatible with the results 
reported in Eq. (14. 3p obtained from the unconstrained fit, which are represented by the 
grey band in the plot. This makes us confident that systematic errors on F coming from 
higher order chiral corrections are smaller than the statistical uncertainties quoted in 
Eq. (|4.3|) . Moreover, it also justifies the topology range that has been chosen at that 
stage to perform the average. 

3 The very weak dependence on the quark mass for v ^ allows us to fix B2 to the expected value 
in the chiral limit. 
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Figure 4: The coefficient B2 of Eq. (|4.2j) as a function of the topological charge \v\. The 
data refer to the smallest quark mass. The results for the lattices el and e3 are slightly 
shifted on the horizontal axis. The dotted line represents the theoretical expectation 
from the chiral effective theory, B2 = 
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Figure 5: The coefficient B\ of the constrained fit Eq. (|4.4p as a function of the topo- 
logical charge \v\. The data refer to the smallest quark mass. The grey bands represent 
the results of Eq. (14.30 . For the lattice e2, results for \v\ = 1 have been omitted due to 
large statistical errors. 
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lat. 


am 


aMp" 


aF^ 


^min / ^ 


d1 


0.02 


209(4) 


0702(8) 


6 




0.03 


0.239(3) 


0.0715(8) 


6 




0.04 


0.267(3) 


0.0727(7) 


6 




0.06 


0.315(3) 


0.0751(7) 


6 


d2 


0.02 


197(3) 


0695(6) 


7 




0.03 


0.2308(24) 


0.0707(7) 


7 




0.04 


0.2604(23) 


0.0719(7) 


7 




0.06 


0.3114(23) 


0.0742(8) 


7 


d3 


0.025 


199(61 


0530(10) 


6 




0.04 


0.242(5) 


0.0551(10) 


6 




0.06 


0.292(5) 


0.0580(10) 


6 




0.08 


0.335(4) 


0.0609(10) 


6 




0.1 


0.375(4) 


0.0637(10) 


6 



Table 4: Results for meson masses and decay constants in lattice units computed in the 
p-regime. t m i n /a indicates the first point in the plateau of aMp eS (t) and aFp eS (t). 



5 Matching lattice QCD with the chiral effective theory: p-regime 



The data obtained for the current correlator C (t) in the p-regime can be compared with 
the expectation of the chiral effective theory at NLO, Eq. Q2.21P ; after symmetrising 
the correlator around t = T/2, we compute the effective pseudoscalar mass aMp eS (t) 
by solving the equation 



(5.1) 



C(t) 


cosh 


\T/2-t)M^ cS (t) 




C(t + 1) 


cosh 


\ T /2-t-l)MV eS (t) 



From aMp cS (t) we then compute an effective pseudoscalar decay constant aFp cS (t) 



FV(t) = 2Z Jt 



C(t)smh[M^ cS (t)T/2 



\|M£ eff (t)cosh [(r/2-i)M£ eff (i) 



(5.2) 



The results for aMp and aFp are obtained from a plateau (for t > t m i n ) and are given 
in Table H 



5.1 Pseudoscalar decay constant 

The quenched chiral effective theory predicts Fp to be volume-independent at NLO 
(see Eq. (|2.22|) ). The volumes of the lattices pi and p2 differ by a factor ~ 3; for each 
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Figure 6: FpTQ as a function of the quark mass (mro/Zs) for the lattices pi, p2, p3. 



quark mass, the results obtained for aFp for the two lattices are compatible within the 
statistical errors 0. 

The quark mass-dependence of aFp is well consistent with a linear behaviour. A chiral 
extrapolation of the form 

aFp = Ci + C 2 {am) (5.3) 
yields for the three lattices the following results: 

Fr = 0.273(4), (pi) 

Fr = 0.272(3), (p2) (5.4) 
Fr = 0.265(6). (p3) 

In Fig. [6] we show the quantity FpVQ as a function of the renormalised quark mass 
mro/Zs and the chiral extrapolation. The chirally extrapolated values are in good 
agreement with each other, indicating that also lattice artifacts are below our statistical 
errors. Alternatively, one can fit Fp linearly in (Mp ) 2 , obtaining the LEC a§ directly 
from the slope: 

Fr = 0.268(4), a 5 = 1.83(13), (pi) 

Fr = 0.269(3), a 5 = 1.69(14), (p2) (5.5) 
4 The correction due to slightly different values of /3 are smaller than the statistical uncertainty. 
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Fr = 0.263(6), a 5 = 1.64(8). (p3) 



The extrapolated values for Fro do not differ significantly from the previous ones. The 
values obtained for 0:5 are consistent with each other and with previous estimates [12,13]. 
The origin of the discrepancy with the results of [10] have been discussed in [13]. 



5.2 Pseudoscalar mass 

We consider the ratio of pseudoscalar meson masses obtained with the lattices pi and 
p2, which have different physical volume (Vx and V2 in the following). This ratio can 
be compared with Eq. (|2.28p . using the topological susceptibility from [58] and Fro = 
0.265(6) from Eq. (pT3|) for the lattice p3, that is 

2 

^0 = ^2 =23.9(2.5). (5.6) 

M P is substituted by Mp 2 , corresponding to the largest volume; once mo and F are 
given, Eq. (I2.28D represents then a parameter-free prediction. The results are reported 
in Fig. [7] (filled circles) for the different quark masses; the empty squares represent the 
parameter-free expectation from the chiral effective theory. The data are compatible 
with the NLO predictions; even though errors are too large to allow for a precise compar- 
ison, NNLO corrections could be significant in this ratio. On the basis of this analysis, 
we can correct our data for finite volume effects according to Eq. (|2.24p and then study 
the quark mass dependence of Mp/(2m). The results are presented in Fig. [HJ where we 
show the data before (top) and after (bottom) the finite volume corrections. Once the 
corrections are applied, M P /(2m) function of the quark mass is well compatible 
with a constant behaviour, without clear evidence for higher order effects. One can use 
anyway the knowledge on NLO chiral effective theory and perform a two-parameter fit 
of the form 

aM p [ w (M p 

{Airy V M 



2m 

with wq and /i given as external inputs. From E\ and E2 we obtain 



2 

+ E 2 M P , (5.7) 



Z( f i) = E 1 F 2 ; 2a s -a 5 = (4vrF) 2 ^. (5.8) 

Using the results for F from the linear chiral fit in m as discussed in the previous section, 
we obtain for the chiral condensate 



l/fj, = 1.5 fm 


S(^)r 3 


= 0.42(4), 


(Pi) 




S(/,)r 3 


= 0.43(3), 


(P2) 




S(/i)r 3 


= 0.41(5), 


(P3) 


I IH = 2.0 fm 


S(/i)r 3 


= 0.47(4), 


(Pi) 




S(/,)r 3 


= 0.48(4), 


(P2) 




S(^)r 3 


= 0.47(6). 


(P3) 



(5.9) 
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Figure 7: Ratio of pseudoscalar masses, (Mp 1 /Mp 2 ) 2 , where V% corresponds to the 
lattice pi and V2 to the lattice p2. The filled circles are the data; the empty squares 
represent the expectation from the chiral effective theory, Eq. (|2.28p . 
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The combination (2as — ot§) can be extracted with very limited accuracy; from this fit 
we obtain errors of order 30-40 %. A more realistic estimation of the errors can be made 
by considering a double-ratio fit 



1 - ^BL log {^-^j + E 2 (M 2 P - M 2 Prcf ), (5.10) 



where one fixes a reference mass m re f and one identifies E2 = 2a% — 05. From this 
one-parameter fit we obtain uncertainties of the order 50-80 % on (2as ~ a b)- It is then 
clear that our precision is too poor to quote a result for this combination. 

6 Comparing results from e- and p- regimes 

After analysing the data in the two regimes, we can now compare the results obtained 
for the leading-order LECs £ and F. 

We first compare the results for Fr^ for the different lattices at our disposal; the data 
from the e-regime, Eq. (|4.3f) . and from the p-regime, Eq. (|5.4f) . are summarised in Fig. [9l 
The agreement is rather good, within one-two standard deviations. The vertical line 
represents the weighted average 

Fr = 0.275(6), (6.1) 

where the error is taken as the largest uncertainty associated to the individual mea- 
surements. We choose this rather conservative error estimate since we know from our 
analysis that systematic effects are below the individual statistical uncertainties. Using 
the phenomenological value ro = 0.5 fm, one obtains F = 108.6(2.4) MeV. 

The results obtained for the quark condensate from Mp/(2m) can be compared 
directly with the ones of a previous work [24] , where the condensate has been computed 
in the e-regime from a finite-size scaling study. The same parameters as for the lattices 
el, e2, e3 have been adopted. The final results quoted in [24] for the renormalisation 
group invariant condensate at the scale L e g = 1.5 fm are 

£ e ff(4rff = 1.5 ftn)rjj = 0.33(3), (el) 

S eff (L efT = 1.5 fin)rg = 0.31(5), (e2) (6.2) 
S eff (L efr = 1.5 fin)rg = 0.29(3). (e3) 

In that analysis it was pointed out that, within the statistical uncertainty, no NLO 
volume-dependence is observed in S e g, as well as no lattice artifacts. We adopted 
Eq. (|2.13|) in order to express all results at the scale 1.5 fm. 

Using the formula in Eq. (|2.30p . we can convert the p-regime results in Eq. (15. 9\\ in 



M 2 P 



2m 
V 2m ref J 



( 
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Figure 8: Quark mass dependence of the squared pseudoscalar mass for lattices pi, p2, 

<M V ) 2 

p3. The plot on the top shows the renormalised quantity 2 ^ r^Zs; on the bottom, 
finite volume corrections according to Eq. (|2.24[) are applied. 
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Figure 9: Summary of the results for Fro obtained in the e- and p- regimes. The 
vertical band corresponds to the weighted average, with error corresponding to the 
largest uncertainty of the individual measurements. 

S e fj(L fj): for the case = L e s = 1.5 fm we get 

E eff (L e ff = 1.5 im)4 = 0.28(3), (pl) 

S cff (L cff = 1.5 fm)r% = 0.285(25), (p2) (6.3) 
S eff (L eff = 1.5 fin)fo = °- 27 ( 4 )- (P 3 ) 

All these results are in good agreement, and are summarised in Fig. [TOl They can be 
converted into the MS- scheme at 2 GeV by using the relation [59] "%jg(2 GeV)/M = 
0.72076. A weighted average gives 

S^(2 GeV) = (292 ± 17 MeV) 3 {L cS = 1.5 fm), (6.4) 

which lies in the same range of previous quenched computations [17-19,23,33,60-63]. 
Also in this case, the error associated to the average is the largest uncertainty coming 
from the single values. 

7 Conclusions 

We have computed the left-left current correlator in quenched lattice QCD, adopting 
two different lattice spacings and two different physical volumes, with quark masses in 
the p- and in the e-regime. We have compared our results with the expectations of the 



21 



L =1.5fm 



0.25 0.3 0.35 



c3 
c2 

el 
P3 
p2 
Pi 



Figure 10: Summary of the results for S e g(L e ff = 1.5 fm)rg obtained in the e- and p- 
regimes. The vertical band corresponds to the weighted average, with error correspond- 
ing to the largest uncertainty of the individual measurements. 



quenched chiral effective theory at NLO. 

In the e-regime we have tested the LO finite-size scaling, which is well reproduced by the 
numerical data. At NLO the effective theory predicts a time and topology dependence 
of the correlator: we have studied both effects and found good qualitative, but a rather 
poor quantitative, agreement with NLO predictions. We should however point out that 
the statistical uncertainty is not good enough to measure the NLO effects with good 
precision. Our results are compatible with the fact that higher order chiral corrections 
could be sizeable at these volumes. We expect these corrections to be more important 
for larger \u\ and indeed the matching seems to work worse there. 
In the p-regime, we have extracted the pseudoscalar mass and decay constant from the 
two-point function. We have checked that the pseudoscalar decay constant at finite 
quark mass Fp is volume- independent, as predicted by the effective theory, and from 
a chiral extrapolation we have extracted the low-energy constants F and a$. We have 
studied the finite- volume effects on the pseudoscalar mass by computing ratios at differ- 
ent volumes. We have observed that after correcting our data for finite- volume effects, 
the behaviour of the quantity Mp/(2m) as a function of the quark mass m is well com- 
patible with a constant one, with no signs of higher order terms, in the range of masses 
that we have considered. Prom Mp j (2m) we have extracted the quark-condensate with 
10% precision. Finally we have compared the results obtained for F from the p- and 
e-regime, finding a good agreement. Moreover, the results obtained for the condensate 
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in the p-regime agree with the ones obtained in the e-regime from a previous study 
based on finite-size scaling. 

This general agreement on the leading-order LECs is a non-trivial test that quenched 
QCD is well reproduced by quenched chiral effective theory at LO. On the other hand, 
a higher statistical accuracy would be needed to test the matching at NLO. 
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A Numerical results 

In this appendix we report the results of the NLO fit according to Eq. (|4.2p of the 
left-left correlators computed in the e-regime. 
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Table 5: Results for the fit in eq. 14.21 for the lattice el, in the time range t/a = 5 — 7 



24 



CLTTl 






Bo 


000316 


o 


08(31 


1(1 4) 




1 


076(9) 

KJ • W 1 KJ y kJ 1 


1.5(6) 




2 


068(3) 

KJ • KJKJKJ \KJ 1 


2.1(6) 




3 


0687(18) 

kj • kj kj i y j- ^ / 


2.7(4) 




4 


0696(19) 

V / • V ' KJ \ i KJ \ -1- ■./ / 


4.0(6) 




5 


0.0706(16) 


3.9(5) 




6 


0.0683(15) 


3.9(5) 




7 


0.0675(20) 


4.5(6) 


000949 





073(20) 

KJ ■ KJ 1 \ W 1 


0.3(1.1) 




1 


075(9) 


1.5(6) 




2 


068(3) 


2.1(6) 




3 


0687(18) 

vy • kj kj kj I y J— ^ / 


2.7(4) 




4 


0696(19) 

V / • KJ KJ \J KJ \ J- '-J j 


4.0(6) 




5 


0.0706(16) 


3.9(5) 




6 


0.0683(15) 


3.9(5) 




7 


0.0675(20) 


4.5(6) 


00253 





067(7) 

KJ * KJ KJ 1 \ 1 / 


0.4(8) 




1 


071(6) 

KJ • KJ 1 _J_ y W I 


1.7(5) 




2 


068(3) 

KJ * KJ KJ KJ \ KJ 1 


2.2(5) 




3 


0687(17) 


2.7(4) 




4 


0696(19) 

V ' • KJ KJ \J KJ \ -L j 


4.0(6) 




5 


0.0706(16) 


3.9(5) 




6 


0.0683(15) 


3.9(5) 




7 


0.0675(20) 


4.5(6) 


00380 





067(5) 

KJ t KJ KJ 1 y ~ J 1 


0.7(7) 




1 


070(4) 

KJ • KJ 1 KJ y 1 


1.9(5) 




2 


069(3) 

KJ . KJ KJ kJ y 'J 1 


2.3(5) 




3 


0686(17) 


2.8(4) 




4 


0696(18) 

V / • KJ KJ \J KJ y -L v.y / 


4.1(6) 




5 


0.0707(16) 


3.9(5) 




6 


0.0680(15) 


3.9(5) 




7 


0.0676(20) 


4.5(6) 


00506 

KJ • KJKJ KJ KJ KJ 


o 


068(4) 

KJ • KJKJKJ \ ~ 1 


1.0(6) 




1 


069(3) 


2.0(5) 




2 


0690(25) 

V ' • V ' KJ \ ' KJ y — " J j 


2.5(5) 




3 


0686(16) 

KJ tKJ KJKJKJ \ J- VV / 


2.9(4) 




4 


0.0696(18) 


4.1(6) 




5 


0.0707(15) 


4.0(5) 




6 


0.0681(15) 


3.9(5) 




7 


0.0676(20) 


4.6(6) 



Table 6: Results for the 



fit in eq. 14.21 for the lattice e2, in the time range t/a = 
25 
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Table 7: Results for the fit in eq. 14.21 for the lattice e3, in the time range t/a = 6 — 10 
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